Bell inequality test with entanglement between an atom and a coherent state in a cavity 
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We study Bell inequality tests with entanglement between a coherent-state field in a cavity and a two-level 
atom. In order to detect the cavity field for such a test, photon on/off measurements and photon number parity 
measurements, respectively, are investigated. When photon on/off measurements are used, at least 50% of detec- 
tion efficiency is required to demonstrate violation of the Bell inequality. Photon number parity measurements 
for the cavity field can be effectively performed using ancillary atoms and an atomic detector, which leads to 
large degrees of Bell violations up to Cirel'son's bound. We also analyze decoherence effects in both field and 
atomic modes and discuss conditions required to perform a Bell inequality test free from the locality loophole. 
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I. INTRODUCTION 

Einstein, Podolsky, and Rosen (EPR) presented an argu- 
ment known as the EPR paradox [1], which triggered the 
debate on quantum mechanics versus local realism. Bell's 
theorem [2] enables one to perform experiments in which 
failure of local realism is demonstrated by the violation of 
Bell's inequality. Various versions of Bell's inequality have 
been developed including Clauser, Home, Shimony and Holt 
(CHSH)'s one [3], and substantial amount of experimental ef- 
forts have been devoted to the successful demonstration of vi- 
olation of Bell's inequality. So far, many experiments have 
been performed to show violation of Bell-type inequalities, 
and most physicists now seem to believe that local realism 
can be violated. 

On the other hand, all the experiments performed to date 
are subject to some loopholes, so that the experimental data 
can still be explained somehow based on a classical (often im- 
pellent) argument. Experiments using optical fields [4—7] typ- 
ically suffer from the "detection loophole" [8], and recent ex- 
periments using atomic states [9, 10] with the maximum sepa- 
ration of ~ 1 m [10], suffer from the "locality loophole" [11]. 
While most of Bell inequality tests have been performed using 
entangled optical fields [4-7], it is an interesting possibility 
to explore Bell inequality tests using atom-field entanglement 
[12-16], particularly for a loophole-free test. In fact, there 
exist theoretical proposals for a loophole-free Bell inequality 
test using hybrid entanglement between atoms and photons 
[14, 17, 18] and relevant experimental efforts [10, 15, 19] have 
been reported. 

In this paper, we study Bell inequality tests with an entan- 
gled state of a two-level atom and a coherent-state field. When 
the amplitude of the coherent state is large enough, such an en- 
tangled state is often called a "Schrodinger cat state" (e.g. in 
Ref. [20]) as an analogy of Schrodinger's paradox where en- 
tanglement between a microscopic atom and a classical object 
is illustrated [21]. Entanglement between atoms and coher- 
ent states has been experimentally demonstrated using cavities 
[22-24]. 

In our study, photon on/off measurements and photon num- 
ber parity measurements, respectively, are employed in order 
to detect the cavity field. We find that when photon on/off 



measurements are used, at least 50% of detection efficiency is 
required to demonstrate violation of the Bell-CHSH inequal- 
ity. One may effectively perform photon number parity mea- 
surements for the cavity field using ancillary probe atoms and 
an atomic detector so that nearly the maximum violation of 
the Bell-CHSH inequality can be achieved. 

The remainder of this paper is organized as follows. In 
Sec. II, we briefly discuss the atom-field entanglement un- 
der consideration and review basic elements of Bell inequality 
tests in our framework. We then investigate the Bell-CHSH 
inequality with photon on/off measurements and parity mea- 
surements, respectively, in Sec. III. Sec. IV is devoted to the 
investigation of the Bell-CHSH inequality test using indirect 
measurements within a 'circular Rydberg atom' -'microwave 
cavity' system. In Sec. V, we analyze decoherence effects in 
both field and atomic modes. This analysis enables us to pro- 
vide quantitative information on the requirements to perform 
a loophole-free Bell test. We conclude with final remarks in 
Sec. VI. 



H. BASIC ELEMENTS FOR BELL INEQUALITY TESTS 

We are interested in testing the Bell-CHSH inequality with 
an atom-field entangled state: 
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W)ac = ^(\e) A \a) c + \§)a\ - <*) c ) , 
V2 



(1) 



where \e)A (I^)a) is the excited (ground) state for the atomic 
mode A, and | + a)c are coherent states of amplitudes +a for 
the field mode C. States (1) for reasonably large values of a 
are considered entanglement between a microscopic system 
and a classical system [20, 25-27]. There have been studies 
on Bell inequality tests with this type of entangled state [20], 
and similar states such as entanglement between an atom and 
a single photon [12-15] and entanglement between coherent 
states [25, 28-36]. Experimental demonstration of state (1) 
has been performed using a system composed of a circular 
Rydberg atom and a microwave cavity field [22-24] . 

In order to test a Bell type inequality, a bipartite entangled 
state should be shared by two separate parties. After shar- 
ing the entangled state, each of the two parties may locally 
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perform appropriate unitary operations and dichotomic mea- 
surements. Violation of the Bell-CHSH inequality can be ob- 
tained by choosing certain values for the parameters of the 
unitary operations. The correlation function is defined as the 
expectation value of the joint measurement 



E(£,0) = (E A (O®E C (]3)\ 



(2) 



where Ea(0 - u a (CFaUa(0 is a dichotomic measurement 
Ta combined with unitary operation Ua(0 parameterized by 
£ , and Ecip) can be defined accordingly. The Bell function S 
is then defined as 

S = \E(£,B) + E{£,B) + E{W) - E{£,B')\ , (3) 

which should obey the inequality forced by local realism, i.e., 
£ <2. The maximum bound for the absolute value of the Bell 
function is 2 V2, known as CireFson's bound [37]. 

An atomic dichotomic measurement can be represented by 
a 2 by 2 matrix 



1 
-1 



(4) 



where we choose the basis as {\e), \g)}. We define the dis- 
placed dichotomic measurement T{() with the atomic dis- 
placement operator D(£) as 



t(0 = D(OrD T (0 



(5) 



with 



D(0 = exp[£d- + -r<5--] = 



cos |£| 



4 sin \^ 



■sin |f| cos 



(6) 
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and < 9 < n and < < 2tt, where & ± are the standard 
ladder operators in the 2-dimensional Hilbert space. We note 
that D{() corresponds to a single qubit rotation for an atomic 
qubit and it can be achieved by applying a Ramsey pulse to 
the atom [38]. We consider measurement F(() for the atomic 
mode A throughout the paper, while some different measure- 
ment schemes are considered for the field mode C. 



III. BELL-CHSH INEQUALITY TESTS WITH 
ATOM-FIELD ENTANGLEMENT 

A. On/off measurement for field mode 

We first investigate the Bell-CHSH inequality with photon 
on/off measurements and the displacement operator for the 
cavity field mode. The displaced on/off measurement for the 
field C is 

00 

OdP) = D C {B){ £ |n><n| - \0){0\)& C (B), (7) 



where £>c(fi) = s\p[Ba c - B*cic\ is the displacement operator 
with the field annihilation (creation) operator ciq (a' c ) and B as 
the displacement parameter for field C. 

We model a photodetector with efficiency r\ by a perfect 
photodetector together with a beam splitter of transmissivity 
sjrj in front of it [39]. The signal field C is mixed with the 
vacuum state |0) v at a beam splitter. The beam splitter opera- 
tor between modes C and v is Bqv = exp[(cos _1 yTTXa^a,, - 
aca\,)l2] [40], where a v (aj) is the field annihilation (creation) 
operator for the ancilla mode v. After passing through the 
beam splitter, the atom-field entangled state | 1 F)ac is changed 
to a mixed state as 



Pac = Tr '' 
" 2 



B C r(mm) AC ®(\o){o\) b ] Cv 

^{\e){e\ ® | y/rja){ y/rja\ + \g)(g\ ® \ - y/rja){- y/rja\ 
e- 2(l -^ 1 \e){g\®\^a){-^a\ 



(8) 



The correlation function with the photon detection efficiency 
r\ is the expectation value of Ta(0 ® @c(fi) f° r state (8) as 



E ((,B; rf) =Tr[pl c f A (0 ® d c (B)\ 



, e -\P\ 2 -\a\ 2 n+2\a\\l3\ V^cosO CQS » 
2 

-> 6 
+ e~ 2|a| cos d> sin - 
2 

- 2 e - 2Ml -W 2+M2 i cos (<p - 2\a\\B\ ^ sin <E>) sin ^, 



(9) 

where a = la-le' ", B = l/^e'*", and O = Op - <D ff with real 
phase parameters <E> ff and Op. The Bell function is immedi- 
ately obtained using Eqs. (3) and (9). 

Using the method of steepest descent [41], we numerically 
find optimized values, |S<9| max , i.e., absolute values of the Bell 
function maximized over variables f, B and B' . We plot 
the results against amplitude lal for various choices of the de- 
tection efficiency from 77 = to rj — 1 (from bottom to top), 
where 77 differs by 0.1 between closest curves in Fig. 1. As- 
suming a real positive value of a, we find that the optimizing 
conditions can also be obtained as 



<r = -, c = 0, b = -b' = 



(10) 



where 1/31 satisfies 



2\B[ 



,2(;;-l)M- 



,-2M|/3| Vfl / 



+ \a\^-e 2M ^(\B\-\a\^. 

(11) 

As expected, the perfect detection efficiency, r\ = 1, gives the 
higher violation up to |Sc)| max ~ 2.61 when \a\ as 0.664. A 
Bell violation of |£o| max « 2.39 (\B \max ~ 2.14) is obtained 
for 77 = 0.8 (77 = 0.6) when \a\ * 0.673 (\a\ * 0.692). 



3 




0.0 0.5 1.0 1.5 2.0 



\a\ 

FIG. 1. (Color online) Numerically optimized values of Bell func- 
tions So with displaced on/off measurements against amplitude a of 
state (1). The detection efficiency ranges in value from rj = (lower 
curve) to rj = 1 (upper curve), with intervals of 0. 1 shown by the fam- 
ily of curves. The horizontal line corresponds to the case of rj = 0.5, 
which coincides with the classical limit of the Bell-CHSH inequality. 



When \a\ = 0, no violation occurs because state (1) con- 
tains no entanglement. As \a\ increases, the Bell violation be- 
comes higher until \a\ ~ 0.7. However, as shown in Fig. 1, 
as \a\ keeps increasing, the degree of the Bell violation de- 
creases towards zero even though the state has larger en- 
tanglement. This result is due to the fact that when \a\ is 
large, the probability of detecting the vacuum for the field 
mode diminishes. Obviously, if photon on/off detection ex- 
cludes one of the two possible results, violation of the Bell- 
CHSH inequality will not occur regardless of the degree of 
entanglement. This is in agreement with a previous result in 
Ref. [30] where the Bell-CHSH inequality with entangled co- 
herent states, \a)\ - a) - 1 - a)\a) (without normalization), was 
considered with on/off detection. 

It should be noted that in Fig. 1, the Bell functions for 
rj = 0.5 overlaps with the horizontal line that indicates the 
classical limit 2. In fact, the photon detector efficiency should 
be higher than 0.5 in order to see a Bell violation as shown 
in Fig. 2(a). Figure 2(b) shows that the optimizing values 
of \a\ are within the range of 0.66 < \a\ < 0.71 for any of 
r\ > 0.5. We also note a previous result [16] that efficiency 
of 0.43 can be tolerated if a different type of Bell inequality 
[42] is used with a nonmaximally entangled state and a perfect 
atomic measurement. 



B. Photon number parity measurement for field mode 

We now consider the displaced photon number parity mea- 
surement for the field mode 

OO 

ftcGS) = £cG6)(Xl l 2 »>< 2 "l - I 2 " + iX 2 " + H^cOS). (12) 
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FIG. 2. (a) Numerically optimized values of Bell function So with 
displaced on/off measurements against detection efficiency ij. The 
local realistic bound, 2, is violated for r\ > 0.5. (b) Plot of optimizing 
values of \a\ with respect to rj. 



Using Eq. (1) and the measurement operators defined above, 
it is straightforward to get 

E n (£,/3)=(t A (0®flc(P)) 

=e- 2 ^~ sin6»cos[4|a|^| sinO - cf>] (13) 
+ e - 2 <M 2 +tf) cos <9sinh[4|a||j6| cos<5], 

and the corresponding Bell function, &n- We present the 
numerically optimized Bell function, |Snlmax, against \a\ in 
Fig. 3, where Bell violation occurs for any nonzero a. Note 
that the atomic displacement operator corresponds to a single- 
qubit rotation for the atomic mode. It was argued that the 
field displacement plays a similar role to approximately rotate 
a coherent-state qubit [30]. If we restrict the atomic displace- 
ment parameters (£ and (') to be real, our test becomes iden- 
tical to the one in Ref. [20] and the result corresponds to the 
dashed curve in Fig. 3. However, it is not sufficient to reveal 
the maximal violation of the atom-field entangled state (1). In 
our numerical analysis, &n is optimized with respect to com- 
plex /3, and f3' that results in the solid curve in Fig. 3. 
Assuming that a is a real positive value, the optimizing con- 
ditions for t>n are found as 

£ = -tt/4, £ = in/4, P = -J3' = i\/3\, (14) 
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FIG. 3. Numerically optimized values of Bell function Sn with dis- 
placed parity measurements against \a\. The solid curve corresponds 
to the absolute values of the Bell function maximized over arbitrary 
f, f, B and B' , while the dashed curve corresponds to those values 
maximized over arbitrary B and B', but real f and f. 
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FIG. 4. (Color online) Schematic of the proposal. The horizontal ar- 
row is to describe the entangled state (1) generation (with R A and C) 
and measurement for atom A. The vertical arrow depicts the indirect 
parity measurement of the cavity field using ancillary atom B. 



where LSI satisfies 



(M-^|)/(M + t0|)=tan4M| 



(15) 



and is nearest to zero. As amplitude \a\ increases, the degree 
of Bell violation rapidly gets larger up to Cirel'son's bound 
2 V2. 



IV. APPROACH USING INDIRECT MEASUREMENT 

In this section, we discuss physical implementations of the 
Bell-CHSH inequality test using displaced parity measure- 
ments in a 'circular Rydberg atom' -'microwave cavity' con- 
figuration. Generation schemes for atom-field entangled states 
(1) have been theoretically studied and experimentally imple- 
mented [38, 43-45]. In the case of a scheme based on the 
off-resonant interaction [38], the required interaction Hamil- 
tonian is 



fi I = h X l(tfa+\)\e){e\-tfa \g){g\\, 



(16) 



and^f = Qr/(46) is the coupling constant determined by the 
vacuum Rabi frequency Q and detuning 6 [38]. As shown in 
Fig. 4, ji /2 Ramsey pulse with phase -n/2 (R A ) is applied to 
a circular Rydberg atom (A) prepared in the excited state \e) A 
[46], which results in an atomic superposition state: \<p-i) A = 
(\e) A - i \g) A ) I v2. Then, a strong dispersive interaction in 
Eq. (16) between atom A and the cavity field produces the 
atom-field entangled state (1) for interaction time t - n/{2x) 
[38]. 

Direct measurements of the light field in the microwave 
cavity are difficult to achieve, while indirect methods for par- 
ity measurements of the cavity field may be more feasible 
[12, 24, 38, 47]. A circular Rydberg atom (B) in Fig. 4 ini- 
tially prepared in state \e) B evolves to a superposition state 
\4>-i>B by 7r /2 Ramsey pulse with phase -n/2 (Rb), and the to- 
tal state is ^Y t ot)ABC = Wac Th e displacement opera- 
tion, D f c (B) = D c (-B), is then applied to the field right before 
atom B enters the cavity, and the same type of interaction as 
Eq. ( 16) between modes B and C follows. One may indirectly 
detect the cavity field by appropriately choosing the interac- 
tion time t = nj(2x) between atom B and the field before 
detecting the atom. The interaction time may be controlled by 
selecting the velocity of atom B. The final measurement for 
atom A, represented by t A (-Tie^ /4), is performed using 7r/2 
Ramsey pulse of phase tt-(/)(Rd) and atomic detector D. The 
measurement on atom B, i.e., Fg(-7r/4), for indirect probing 
is performed with the help of n/2 Ramsey pulse with n phase 
(Rd') and atomic detector D' . The measurement operator is 
then represented as 



f BtC (j3, t) = <U B ,c(p, t) O b , c <U b .c(B, t), (17) 



where 



B ,C = l\+){ + \-\-){-\]B®tc, 

and |±> = (|e> + \g))l V2. The correlation function E(£,8, t) 
\Ta(0 ® t B>c 03, o) is calculated using state | v P tof )ABc as 
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p, t) =^cos0 e (w2+ ^- 2| '* |cosO,( - 1+cos2 ^cos[(k| 2 + 



2\a\\B\ cosO) sin2^f] 



_ CQS 0e (l«| 2 +l8| 2 +2|a|l8|cosO)(-l+cos2^) 

2 



cos[(|a| 2 + |y0| 2 + 2|or||/?| cos 0>) sin2^f] 



H — sin (fe 
2 

1 

H — sin 8e 
2 



W : -^-(N 2 -^| 2 )cos2^-2W Ws inO S in2^ cos[ ^ _ ( | ff |2 _ \ B \2 )sin 2 X t + 2\a\\p\ sinO(-l + COS 2#t)] 

-|«| 2 -^| 2 -(N 2 -^| 2 )cos2^ + 2M Ws in* s i„2^ cos[ ^ + ( | a |2 _ p2 )sin2; ^ + 2 |a|^| sinO(-l + COS 2^f)] 

I 



(18) 
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Cylindrical cavity 
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T> C (J3)' applied to cavity field C 
when atom A passes here. 

FIG. 6. (Color online) Sideview of the atom A's path with intervals of 
time. Each interval denotes an amount of time required for atom A to 
pass through the region related with atomic velocity v. We note that 
the distance I = v X {U + h), which corresponds to the length of the 
cylindrical cavity, is a crucial factor in a loohole-free Bell inequality 
test. 



FIG. 5. Numerically optimized values of Bell function St with 
indirect measurements against |or|. The result is found to be identical 
to the one using direct parity measurements shown as the solid curve 
in Fig. 3. 



tative requirements to perform a loophole-free Bell test. 



A. Decoherence effects in the cavity-atom system 



and the Bell function, St, is accordingly obtained. As ex- 
pected, the optimizing conditions for |Srlmax are identical to 
those for |S<9| max in Eqs. (14) and (15) with an additional con- 
dition, t — t' — n/2x- Our numerical study confirms that 
the optimized Bell function | St I max plotted with the above- 
mentioned conditions in Fig. 5 exactly overlaps with the solid 
curve in Fig. 3 as shown . This result is due to the fact that 
the indirect measurement (17) is basically equivalent to the 
displaced parity measurement (12) on the cavity field when t 
is chosen to be n/(2x) [47]. I.e., the measurement on atom B 
in the basis (|+) , |-)} after the interaction time t = n/(2x) is 
equivalent to the parity measurement on the cavity-field. In 
fact, it can be shown that the correlation functions (18) with 
t = nj{2x) an( l (13) are identical. Of course, if we restrict 
(, to be real, the optimized plot of the Bell function ISrlmax 
approaches the dashed curve in Fig. 3. 



V. DECOHERENCE AND LOOPHOLES 

It is not difficult to predict that decoherence effects due to 
the cavity-field dissipation and the spontaneous emission of 
the atoms will obstruct Bell violations. This is particularly 
important when one intends to demonstrate a Bell violation 
free from the loopholes. In this section, we consider deco- 
herence effects with realistic conditions for the Bell-CHSH 
inequality test using parity measurements and suggest quanti- 



There are two main effects that cause decoherence in our 
Bell inequality test, i.e., spontaneous emissions from atoms 
and cavity field dissipations. In the atom-cavity system under 
consideration, one (or both) of these two effects may occur. 
The master equation which determines the time-evolution of 
the density operator, p{t), under the atom-field interaction with 
spontaneous emissions and cavity dissipations is 



dp(t) 1 - 
at in 

with the Linblad decohering term X. defined as 

Hp =K(2apa^ — a^cip — pa'' a) 

+ y{2&-pcr + - & + &-p - p& + &S), 



(19) 



(20) 



where k is the dissipation rate of cavity field, and y is the 
spontaneous emission rate. 

It is known that the spontaneous emission rate of an atom 
can be significantly reduced by engineering the shape of the 
cavity that contains the atom [48, 49]. A complete inhibition 
of spontaneous emission was suggested using a cylindrical 
metal cavity with a diameter shorter than 1.8412e/wo, where 
a>o is the transition rate between atomic states \e) and \g) and 
c is the speed of light [49]. For our setup, the transition rate 
can be taken from Ref. [24] as ioq = 5 1.1 GHz. This means 
that the diameter should be smaller than 3.44mm that is ex- 
perimentally achievable. As seen in Fig. 6, a long cylindrical 
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FIG. 7. (Color online) A timeline for decoherence with dynamical 
parameters related in each regions (from left to right). The top line is 
for atom A, the middle for cavity field C, and the bottom for atom B. 
The times when Ramsey pulses are applied are described as vertical 
dashed lines. We consider a Ramsey pulse application as an instant 
event as a Ramsey pulse lasts as short as 1 jis order [45]. Regions 
are differently hatched depending on the types of dynamics. In the 
diagonally hatched regions, atoms A and B travel in free spaces with 
the spontaneous emission rate y before and after Ramsey pulses as 
shown in Fig. 4. In the cross-hatched region, atom A travels in a 
cylindrical cavity with the inhibited spontaneous emission rate y' . In 
the vertically hatched region, the cavity dissipation with rate k occurs 
in the cavity (C) field. The horizontally hatched regions correspond 
to the dynamics of the atom-field interaction Hj in the main cavity 
C together with spontaneous emission y c and cavity dissipation k. 
Abbreviations C, D, D', R A , R B , R D , and R D i are consistent with 
those in Fig. 4. 



cavity may be used between cavity C and Ramsey zone Ro to 
inhibit spontaneous emission. 

The spontaneous emission rate y c inside the cavity C in 
Fig. 4 is also generally different from the spontaneous emis- 
sion rate yo in the vacuum. It is known that y c can be calcu- 
lated by approximating the cavity in the one dimension while 
considering the effect of the atomic motion as described in 
Ref. [50]. In our case, y c = 4.08 Hz is obtained based on the 
result of Ref. [50] from the spontaneous emission rate in the 
vacuum, yo = l/(27o) (To = 36 ms is the atomic life time in 
the vacuum [38]) and related realistic parameters in a recent 
experiment [24] . 

Considering the discussions above, we present a timeline of 
decoherence effects in Fig. 7 together with time intervals re- 
quired to pass through certain parts of the apparatus as follows 
(also depicted in Fig. 6): t\ is a half of the time required for 
an atom to pass through a cavity used for Ramsey pulse ap- 
plication, ?2 is a half of the time required for an atom to pass 
through Ramsey pulse and the main cavity (C) without cav- 
ity waist, ?3 for an atom to pass through the main cavity(C)'s 
waist (n/2x), U for atom A to pass through the long cylindri- 
cal cavity before the field displacement operation on the cavity 
field, ts for atom A to pass through the remainder of the long 
cavity after the field displacement operation, and ?6 for atomic 
detection at D or D' . 

Let us first consider the pathway of atom A, which corre- 
sponds to the top line of Fig. 7. Atom A undergoes sponta- 
neous emission before and after the Ramsey pulse Ra with rate 
yo (diagonally hatched part). Atom A then interacts with the 
cavity field with dissipation rate k under spontaneous emission 



(jc), which is represented by the horizontally hatched part. 
After the atom-field interaction, atom A passes through the 
cylindrical cavity experiencing inhibited spontaneous emis- 
sion {y'). Finally, atom A comes out of Ramsey pulse Ro 
experiencing spontaneous emission (yo), and is registered at 
detector D. In the mean while, cavity field C which have inter- 
acted with atom A undergoes field dissipation (k) while atom 
A is passing through cylindrical cavity. Then, cavity field C 
begins to interacts with atom B under spontaneous emission 
(y c ) and field dissipation (k) after displacement operation on 
it. Atom B, used for an indirect measurement, experiences 
spontaneous emission (yo) around Ramsey pulse Rb, interac- 
tion with the cavity field (C) with spontaneous emission (y f ), 
and spontaneous emission (yo) before detection D'. 

Here, we take the photon storage time Tq = 0.13 s (k = 
l/(27c)), Q = 2tt-49 kHz and 6 = 2tt-65 kHz (% = £l 2 /(46) * 
58kHz) from recent experiments [24]. The solution of the 
master equation for the cavity dissipation alone with Hi, was 
examined in Ref. [51]. In Appendix, we obtain the solution of 
Eq. (19) and find an explicit form of the density operator and 
the correlation function. The Bell function can be constructed 
using the correlation function in Eq. (B13) of Appendix. Note 
that we have assumed perfect Ramsey pulses during the pro- 
cedures. Considering cavity dissipation, we employ the same 
optimizing conditions (14) except that \B\ is chosen to be the 
values that satisfy 



and is nearest to zero. 



tan(4|ff|e 



-K(ti+h) 



(21) 



B. Bell violation and separations under practical conditions 
without a cylindrical cavity 

Let us first consider Bell violation depending on the sep- 
aration / = v x (?4 + ts) between both parties without using 
a cylindrical cavity (thus y' = yo). We choose some practi- 
cal time-interval parameters as t\ = 80.0 /ys, tz - 166.5 ps, 
?3 = 27.1 yus, ?6 = 20 fis and velocity of an atom v = 250 m/s 
[24, 52]. The Bell function with several choices of / are plot- 
ted in Fig. 8. The Bell function approaches the value near 2.7 
when / = 0.1 (meter), but it decreases as / gets larger. Clear 
Bell violations appear for / < 2(meter), however, this is in- 
sufficient for a space-like separation as we shall discuss in the 
next subsection. 



C. Requirements for a Bell test free from the locality loophole 
with a cylindrical cavity 

In principle, a Bell test free from the locality loophole can 
be performed using a long cylindrical cavity with a low spon- 
taneous emission rate (y') and the main cavity with a low dis- 
sipation rate {k). In order to close the locality loophole, the 
measurement event for atom A should not affect the measure- 
ment event for the cavity field C, and vice versa [11]. In other 
words, the measurement event for atom A should be outside of 
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FIG. 8. The Bell function under realistic conditions discussed in 
Sec. VB are plotted with optimizing conditions in Eqs. (14) and (21) 
for several different cases of separation /. As the separation / gets 
larger, the maximum values of the Bell function decrease. The deco- 
herence effects become heavier as \a\ gets larger. 



the "back light cone" from the detection event D' in Fig. 4. In 
the same manner, the measurement event for the cavity field C 
should not be in the back light cone from the detection event 
D. For simplicity, let us first suppose that each measurement 
process takes place at a single location (D and D'). In our 
Bell test, the time t& required to measure atom A is smaller 
than the time required to measure field C (tc) due to the in- 
direct measurement scheme for field C. We assume that the 
measurement event for the field C precedes to the measure- 
ment event for atom A by T (the opposite case will require a 
longer separation between the two parties). Then the condi- 
tions required to close the locality loophole are 



d > c(T + t A ), 
d > c(t c - T), 



(22) 



where d is the distance between D and D' and c is the speed 
of light. 

In order to apply the locality-loophole-free conditions (22) 
to our Bell test setup in a more rigorous manner, one needs 
to consider locations of the local measurement elements. In 
Fig. 7, one can find that the measurement time for atom A (t A ) 
consists of the times for Ro(t\) and D (t(>) and that for the field 
(?c) consists of the times for C fe), R' D + ft), and D' (t^). A 
measurement event for each party actually does not take place 
at a single location, and both of the measurements are not even 
on a straight line. Therefore the distance d in Eqs. (22) needs 
to be replaced with the distances from the final detector of 
one party to the location where the measurement of the other 
party begins. A careful consideration leads to the conclusion 
that the following inequalities should be satisfied: 

+ t 4 + ts + t\ + t 6 ) > c(t 5 + t\ + f 6 ), 

v Vfe/2 + t 2 + h+ t 6 ) 2 + (f 3 /2 + U + t 5 ) 2 (23) 
> c(h + t 2 + h + t 6 - t 5 ). 

Using the feasible values of t\, t 2 , ?3, t(, and v in the previous 
subsection, we find the minimum values t 4 = 236.0 s and t$ - 
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FIG. 9. (Color online) Contour plots of the Bell function with respect 
to photon storage time Tc in the main cavity and amplitude \a\ of the 
entangled state. The atomic life time in cylindrical cavity r atom is 
fixed at 1000, 2000, 4000 and 00 (seconds). The minimum distance 
condition / = 52.99 (km) for a loophole-free Bell test was assumed. 
Higher inhibition of spontaneous emission in the cylindrical cavity 
reduces the required photon storage time in the main cavity. 



96.8 //s with which the equalities hold for Eqs. (23). Then, 
the minimum distance required for a Bell test free from the 
locality loophole is found to be / = 52.99 km [24]. 

We finally consider conditions of the atomic life time 
T'atom = l/(2y') and the photon storage time T c = 1/(2a:) re- 
quired for a Bell test free from the locality loophole. In Fig. 9, 
we plot the Bell function constructed using Eq. (B13) in Ap- 
pendix with respect to the photon storage time in the main 
cavity and amplitude \a\ of the atom-field entanglement. Here, 
the extended lifetime of the atom in the cylindrical cavity was 
assumed to be r atom =1000, 2000, 4000, and oo (seconds). 
The distance I was assumed to be the minimum distance re- 
quired for a loophole-free Bell test (52.99 km). For exam- 
ple, when r atom = 2000 (seconds), the photon storage time 
Tc ~ 1 160 (seconds) at |ff| ~ 0.47 is required to see a Bell vi- 
olation. If complete inhibition of the spontaneous emission in 
the cylindrical cavity is possible, (i.e., T atom = oo), Tc ~ 590 
at \a\ ~ 0.3 is required. Obviously, the stronger inhibition of 
the spontaneous emission in the cylindrical cavity relaxes the 
requirement of the photon storage time in the main cavity to 
see Bell violations. However, it still requires at least a few 
hundreds of seconds for the photon storage time to demon- 
strate a loophole-free Bell violation, while it is only about 
0.13 s at present [53]. It would also be extremely challeng- 
ing to build a long cylindrical cavity that strongly inhibits the 
spontaneous emission of atom A during such a long life time. 
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VI. REMARKS 



1. Spontaneous emission for atom 



We have investigated Bell-CHSH inequality tests with en- 
tanglement between a two-level atom and a coherent-state 
field in a cavity. In order to detect the cavity field for these 
tests, photon on/off measurements and photon number parity 
measurements, respectively, have been attempted. When pho- 
ton on/off measurements with the perfect efficiency are used, 
the maximum value of the Bell violation is So ~ 2.61 at 
\a\ ~ 0.664. In order to see a violation of the Bell-CHSH 
inequality, at least 50% of detection efficiency is required. 
When photon parity measurements are used, the value of the 
Bell-CHSH violation rapidly increases as a gets larger, and 
it approaches CireFson's bound for a » 1. Although pre- 
cise direct measurements of cavity fields are experimentally 
difficult, photon number parity measurements for the cavity 
field can be effectively performed using ancillary probe atoms 
and atomic detectors. We have fully analyzed decoherence 
effects in both field and atomic modes and discuss conditions 
required to perform a Bell inequality test free from the locality 
loophole. 

Our proposal may be considered an attempt to analyze 
a Bell inequality test using entanglement between a micro- 
scopic system and a mesoscopic classical system. Since 
atomic detectors are known to be highly efficient [54], it may 
also be a reasonable target to perform this type of experiment 
in a way free from the detection loophole. In principle, a Bell 
inequality test free from the locality loophole in our frame- 
work using atom-field entanglement may be performed using 
a long cylindrical cavity for the atom with a low spontaneous 
emission rate [49]. However, our analysis shows that it would 
be extremely demanding to perform a Bell inequality test free 
from both the locality and detection loopholes in this frame- 
work since the main cavity for field with a low dissipation rate 
would be necessary together with a long cylindrical cavity. 
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Appendix A: Solutions of the Master Equation for Matrix 
Elements 



A density operator of a two-level atom, p A (t), can be ex- 
pressed as a matrix form 

where pA,ij(f) = {i\pA(f)\j). When an atom with a initial 
density matrix, Pa(0), goes through the spontaneous emission 
process for time f, its density matrix is straightforwardly ob- 
tained using Eq. (19) with^ = and k — as 



p A (t) =s A (r,t)[pA(0)] 



e- yt PA, eg (0) 



(A2) 



e 7 'p A ,ge(0) PA,gg(0) - PA,ee(0)(e 



Zyt _ 



1) 



where superoperator S(y, f) is defined for later use. 



2. Dissipation for cavity field 

In order to find the time evolution of the coherent-state part 
the density operator, it is sufficient to find the time evolution of 
an operator component \fx) (v|, where \p) and |v) are coherent 
states of amplitudes fi and v. This solution for time t under the 
master equation (19) with x — and y = is well known as 
[55, 56] 



exp[- 



(M 2 + H 2 -2v^)(l-exp(-2^)) i 



hie 



'}<• 



(A3) 



3. Atom-Held interaction with spontaneous emission and 
cavity dissipation 

The density matrix p(t) for an atom-field state can be con- 
sidered in a 2 x oo dimensional space, since we assume a 
two-level atom. It is possible to decompose the master equa- 
tion (19) in \\e) , \g)} basis with the density matrix elements 
Pdj(t) = (i\p(t) [/)■ We then obtain equations 



dt 



PC.ee =LeePC,ee ~ 2yp c , £ 



ix[a l a,p c ,ee] + K{2ap aee a' - a' ap c , e 

2yPC,ee, 



■ pc, ee CVa) 

(A4) 



We first find general solutions of the master equation (19) 
for three types of decoherence processes step by step, i.e., 
spontaneous emission of an atom, cavity dissipation, and 
atom-field interaction with spontaneous emission and cavity 
dissipation. 



~rf~ t PCgg -£ggPC,gg + 2yPC,ee 

=ix\ja*a,pc r gg\ + K(2apQ gg a' - a^apc, gg - pc,gga T a) 

+ tyPCec 

(A5) 



9 



d_ 

dt 



PC.eg =LegPC,eg ~ IXPC.eg ~ JPC,eg 

= - ixia'a + \)pc,eg ~ iXPc,e g a*a 
+ K{2ap c ,egCi ] - a f ap c ,eg -pc,e g a ] a) - yp c , e 



(A6) 



■J t PCge =£gePC,ge + IXPCge ~ JPCge 

=iXPc, g e(a ] 'a + 1) + ixa f ap c , g e 

+ K(2ap c ,g e a* - a'ap Cge - p Cge tf&) - yp c , ge - 

(A7) 

We define the following superoperators for simplicity: M — 
a'a ■ , P — ■ a^a , and J' - a ■ fi' . Then £ ee , £, gg , £, eg , 
and X g e can be expressed as, 



L ee = 2kJ -rM- r*P, 
L eg = 2kJ -rM-rP, 



(A8a) 
(A8b) 



where r = K + ix, and £. gg and £, ge are obtained by substituting 
X with ~x in £^ e , and £, eg , respectively. A master equation 
of the form dp/dt = £p + cp, where c is a constant and _£ is 
a superoperator, can be solved with a usual exponential form 
exp [(£ + c)t]p. 

a. Solution for pc,ee 
The solution of Eq. (A4) is 

p C ,ee(t) = explilee ~ 2y)t]p Qee (0) = e^' e 1 "' p C ,ee(0) . (A9) 

where the factorization can be done by the similarity transfor- 
mation [57]. Now we need to factorize er- at . This is solved 
with an ansatz (a technique can be found in Ref. [56]) 

pc.eeit) = exp[-2yf] exp[(-rM-r*P)t] exp[/(/)2*JTp c , ee (0), 

(A10) 

where f(t) — (1 - e~ 2l ")/(2i<). For an initial state pcw(O) = 
^>(v|, 

Pc,eM = exp[-2yf + &(k, 0,ju, v, t)] ^ie~ rt ) (ve" rt | (Al 1) 
with 

®(K,x,»,v,t) := -^(|v| 2 + \p\ 2 )(l - e- 2 "') + -(1 - e - 2 '» 
2 r 

(A12) 



b. Solution for pc,gg 

In order to solve Eq. (A5), we first assume y — 0. A homo- 
geneous solution is obtained from Eq. (Al 1) by substituting^ 
with - x as 



Then it is obvious to see that the general solution pc, gg (t) with 
y + is 

Pc,g g (t)=p h Qgg (t) + 2y f dt'p Q ee(t') 

= exp[0(K, 0,//, v, f)] ^e" r * f ) (ve" r *'| 

+ 2y J~ <if' exp[-2yf' + 0(k, 0,/i, v, f')] (ve~"'| . 

(A14) 

c. Solution for pc,ge 

The solution of Eq. (A7) is 

pc, ge (t) = exp[(lge + ix~ y)t]pc, ge (0) = e^' e l *>' p c , ge (0) . 

(A15) 

Factoring er-* 1 with an ansatz 

p c ,ge(t) = exp[(ix-y)t] exp[(-rM-r*P)t] exp[g(t)2i<J]p Qee (0), 

(A16) 

where g(t) = (1 - e - 2r * f )/(2r*). Forp c ,^(0) = |/i> <v|, 



Pc,ge(t) = exp[(i X - y)t + ®(k, -x,H, v, t)] \pe ''*') (ye r '\ 



(A17) 



d. Solution for pc.eg 

The solution of Eq. (A6) is obtained from Eq. (A17) by 
substituting x with ~x as 



Pc, eg (t) = exp[(-i* - y)t + ®(k,x,P, v, f)] \pe ''') (ye r ''\ 



(A18) 



Appendix B: Derivation of the density matrices for atom-field 
entanglement and the correlation function 

1. Atom-field entanglement generated under decoherence 
effects 



First, atom A initially prepared in \e) A undergoes the spon- 
taneous emission for the time t \ . After applying the first Ram- 
sey pulse, R A = D A (-i7r/4), explained in Sec. IV, atom A 
again undergoes the spontaneous emission for time t2- Using 
Eqs. (A2) again, it becomes 



S A (yo, t 2 )[D A (-in/4-)\{s A (y , h)[{\e){e\) A ^ A {-i7tlA)] 

i e -2y f2 (_| + ie -2r^ e -yx>h \ (Bl) 

a - f e - 2 vo")e-v«'- a - l -e- 2 y° h -) 



Then atom A interacts with the cavity field C prepared in state 
\ia) c . Using Eqs. (All), (A14), (A17) and (A18), we find 
the state after the interaction time t$ as py c with its matrix 
elements: 



p h Cgg (t) = exp[0(*, 0,fx, v, t)] )fie-") (> 



(A13) 



PAC.e 



l e -2yoh~2y cti U 



(B2a) 
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P ( AC,e g = (-5 + '^ 2ro ")exp[- ro f2 + (-«■ - 7^3 
+ ®(k,x, a, a, f 3 )] \iae~ r '^ (iae~ r ' h \ , 



P?c, g e = (j - '^ 2ro ' 1 ) exp[-rof 2 + (fr - r,)f3 
+ 0(k, a, a, f 3 )] |j'ae~ r *' 3 Wjae~'' t3 | , 



(B2b) 



(B2c) 



(B2d) 



2. Atom-field entanglement after traveling for the spacelike 
separation 

We now derive the total density matrix right before D f c (B) is 

applied. The state p A ^ undergoes spontaneous emission inside 
the cylindrical cavity and dissipation inside the cavity field C. 
The calculation can be done using the results in Sec. A 3 with 
X — 0. Then the state becomes p A l, where 



B, which is in the same state as that of atom A in Eq. (Bl), 
goes into the cavity field of state ff AC . The atom-field interac- 
tion Hi with the coupling constant x occurs between atom B 
and field C for time ? 3 . When solving the master equation, it is 
convenient if one notes that the field-part of state ff AC can be 
expressed by coherent-state dyadics such as \fx)(v\. If the com- 
ponent of the cavity field, initially prepared as |/z)(v|, interacts 
with an atomic state (Bl) for time f 3 , the resulting density op- 
erator element is obtained as Cl B c(M> v > h) with 



&BC,ee(H, V, f 3 ) = ~e 



(B4a) 



&BC,eg(M, v,h) = {-- + ie 2yo,] ) exp[-y f 2 + Hx ~ 7c)h 



+ ®f>,*,/i,v,* 3 )] ]fie-^)(ve- r ^ 



(B4b) 



&BC,ge(M, v, h) = G - ie 2ro " ) exp[-y f 2 + (ix ~ 7c)?3 



®(K,-x,fx,v,t 3 )]\fxe- r% ){ve-^\ 



(B4c) 



pfcee = -e- 2yah - lyA - ly,ti \iae-**-*) (/ae^' 3 ^' 4 | , (B3a) 



P ( AC,e g = (~2 + ie ~ Z70tl ) ex Pt-ro?2 + (~ix ~ Jc)h - y'h 
+ ®(k,x, a, a, f 3 ) + ©0, 0, ae- r '\ ae~ r " h , h)] 
\iae- r ' 3 - Kti )(iae- r ' t3 - Kli \, 



i 



(B3b) 



pfcge = ( 2 - ie ~ 2rotl ) ex P[-r'2 + dx - jc)h - y'h 

+ ®(k, -x, a, a,t 3 ) + ®(k, 0, ae~ r '' 3 , ae"''' 3 , f 4 )] (B3c) 
\iae- r% - Kk )(iae- rt >-*' 4 \, 

PAc? g = (1 - \e-^)\iae-' l >-«)(iae-''>>-«\ 



+ 2y 



' C dt-e 
Jo 2 



-2y oh -2y ch -2y'< \ iae -rt 3 -«t} (i ae ~rt 3 - 



(B3d) 



fi*toGi.v,f 3 ) = (1 - ie- 2 ^)e & ^"^ \pe-^)(ve-''' 



2y, 



f 3 j 1 

Jo 2 



-2y o '2-2y r r+0(*,O,/j,v,O 



|jue- rt )(ve- rt | 



(B4d) 



We used Eqs. (All), (A14), (A17) and (A18) again to find 
Eqs. (B4a)-(B4d). Therefore, p{ c interacts with atom B and 
evolves to 



Pabc =(\e)(e\)A ® pBC.ee + (\e)(g\h ®Psc, t 



(\g)(e\)A ®PBC,ge + (\g)(g\)A ®PBC,gg, 



where 



pBC,ee ——e 



-2y 1 2 -2y c (3 -2/ (I4 +f 3 ) 



(B5) 



(B6a) 



Clucdae ~" 3 ~ ai -BJae-" 3 -"' 4 -B,t 3 ), 



PBQeg =(-- + ie 2yoh ) exp[-y f 2 + (-ix ~ Jc)h - y'(h + f 3 ) 
+ ®(k,x, a, a, f 3 ) + ®(k, 0, ae~"\ are~ r *' 3 , t A ) 
- 2ie- K( ' 1+,3) sinCrf 3 )Im(or*j8)] 
ClBcdae-^-^ -B,iae- r ""- K 'i -B, f 3 ) 

(B6b) 



and here the subscripts are consistent with the previous ones. 



3. Effects with atom B for indirect measurement 



After applying the displacement operation D C (B), the total 
state becomes ff AC = T)£(B)p AC l)c(j3). Now, the probe atom 



PBQge =(~ - ie 270,1 ) exp[-y f2 + (ix ~ 7c)h - Y(h + h) 
+ ®(k, -x, a, a, / 3 ) + ®(k, 0, ae - ''*' 3 , ae - " 3 , f 4 ) 
+ 2/e" ,f( '- 1+r3) sin(^f 3 )Im(a , ' J 6)] 
n B c(iae- r ' l3 - Kk -p,iae- rt *-" h -B, f 3 ) 

(B6c) 
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PBC.gg = (1 - -e~ 2yah )Q. B c{iae~ r% ~' <ti -8, iae~ r% ~ Kti -J3, t 3 ) and operator U B (jU, v, t 3 ) = TT C Cl BC ifi, v, t 3 ) is determined as 



2y c ^ d^-e-^-^Cl 

r' 1 1 

2/ l *2 



BC {iae- r, - K,i -p,iae-"- K,i - B,t 3 ) 



c ,-2-y t 2 -2y c t3-2y't i 



Q. BC (iae~ 



■ B, iae 



X(l -e- 2y,i ) 



(B6d) 



Spontaneous emissions of atom A and B that may occur after 
this point shall be taken into account when we derive the cor- 
relation function. Since the field state is not considered any 
more from this point before the final measurements, the cav- 
ity dissipation can be ignored. By tracing out the cavity field, 
we get 



Pab = TtcPabc = ^ (\i)(j\)A ® o-Bjj, 



(B7) 



where 



0-B.ee =Tr C p BC 



1 



, -2y Q h -2y c r 3 -2y ' (r 4 +f 3 ) 



(B8a) 



\5 B (iae 



-rti -Kit 



a- B , eg =(~2 + ie ) exp[-y f 2 + Hx ~ Jc)h - j'(h + t 3 ) 
+ &(k,x, a, a, t 3 ) + ®(k, 0, ae~ n% , ae~'*' 3 , f 4 ) 
- 2ie- K( '* +li) sin0t-f 3 )Im(a*j8)] 
■U B (iae- r '^ - B, iae~ 1 * h ~ Kti - B, t 3 ), 

(B8b) 



d- B , g e =(~ - ie 2rofl )exp[-7 f2 + (ix - J c )h - y'(k + t 3 ) 
+ ®(k, -x, a, a, t 3 ) + ®(k, 0, ae~ r ' h , ae~ rh , f 4 ) 
+ 2ie- K{u+h) sm(xt 3 )lm(a*B)] 
, & B (iae- r% - Kti -B,iae-"i- ai -B,t 3 ), 



(B8c) 



d-B.gg = (1 - -e 



-*ro>2-2rctjj Biiae -rt-Kt A 



l5 B ee(H V t 3 ) = _ e - 2 yo'2-2y r f3+0fcO,/i,v,r 3 )-i(W 2 +|v| 2 -2/jv*)exp(-2<f/ 3 ) 

(B9a) 



U B , eg (M, v, h) = (- 2 + ie 2yoh ) exp[-y f 2 + Hx ~ Jc)h 



+ ®{K,x,n,v,h)- -( 



|v| 2 - 2/zvV 2! *' 3 )<r 2 *' 3 ], 



UB, ge (v, v, t 3 ) = (- - /e - yoh ) exp[-y f 2 + (fr - Jc)h 



+ ®(K,-X,H,V,t 3 ) 



+ \v\ 2 -2iiv*e 2ixh )e- 2Kh } 



(B9b) 



(B9c) 



U B , gg (M,y,t3) = (1 - I e -2r'2 )e 0(^.v, f3 )4( W 2 + |vP-2^*)exp(-2,r 3 ) 



+ 2y, 



(B9d) 



4. Decoherence right before final measurements and the 
correlation function 



We now consider the last measurement process for both 
parties. Atom A experiences spontaneous emission for time 
?5 - t 3 with rate y', then atomic displacement operation 
Dj^—e^nlA) is applied. After the displacement operation, 
atom A evolves again under the spontaneous emission for time 
t\ with rate yo. We define superoperator X to describe this 
process as 

XA(y,yoJ5-t 3 ,h,cf>)[pA] 

= <S A (y ,ri)[^(- e -'V4){>S A (r',f5 - t 3 )\p A ]}A A (re-»n/4)] 

(BIO) 

Atom B undergoes spontaneous emission for time f 2 with rate 
yo, and displacement operation D B (-n/4) is applied. Then, it 
experiences spontaneous emission for time t \ with rate yo just 
before the final measurement. This process can be expressed 
as 



X B (yo,yo,t2,ti,0)lp B ] 

= S B (y , h)[6l(-n/4){3 B (y , t 2 )[p B ]}t>B(-n/4)] 



(Bll) 



2y c f 3 dt\e- 2y ° h -- 2y < 'i3 B (iae- r, - K,i - B, iae'^ - B, t 3 ) 
Jo 2 

2y' f 4 dt-e- 2yah -- 2yA - 2y ''i3 B (iae- rh - K ' - B, iae-"'-" - B, t 3 ) 
Jo 2 

-BJae-'"-" 1 -B,t 3 ) 

2y ' h ), 



The final density operator used to obtain the correlation func- 
tion is then obtained using state p^ B in Eq. (B7) with Xa and 
X B as 



+ _ e -2rotz-2 7c t 3 -2y k Q B ( 



x(l 



Pab = X A {y',y(hh - t 3 ,tu4>)®X B {y ih y Q J 2 ,tuQ)\pAB\- 

(B12) 

The correlation function is obtained as the expectation value 
of dichotomic measurements (4) performed by both the par- 
(B8d) ties: 
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where 



EQfi.fi, t u t 2 , t 3 , k, t 5 ) = Tr[p« n /f A 9 t B ] 

= 1 -2y , 2 -2y ct} -2y>, 4(e -2y n _ _ g -2/ f3 + ^V*^^, A lf t 3 ) + £ &(Al, A 2 , f 3 ) + Z*&(A 2 , A lt f 3 ) 



2 



+ (e -2w, _ _ L- 2 »)^(A 2 ,A 2 ,? 3 ) + 2 rc V dte-^-^^Qae-"-^ -fijae'^ -fi,t 3 ) 

Jo 2 

+ 2y' f 4 dt-e-^- 3 *° t >- 2 *' , &,(iae- ,,s -" - fi, iae^-* 1 - fi, f 3 )}, 
Jo 2 

v, f 3 ) = {^B, ee Qi, v, h) + V B , gg (M, v, f 3 ))( e - 2y °" - 1) + CU B ,e g (M, v, h) + U^Ou, v, t 3 ))e-^- 2 ^ , 

+ i e -2yuh^ e -yo(t2+2n)M-ix-yc)'3-y'('4+t5)+®(K,X,a,a,t } )+@(KX)^^ 



(B13) 



(B14) 



I 

Ai = iae~ rt3 ~*' 4 - [5 and A 2 = iae~ r '' 3 ~ 1 " 4 - fi. Using this cor- 
relation function, one can eventually construct the Bell func- 
tion using Eq. (3). 
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